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S_i 1 ABSTRACT 



The spectroscopic Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7) galaxy sam- 
ple represents the final set of galaxies observed using the original SDSS target selection cri- 
teria. We analyse the clustering of galaxies within this sample, including both the Luminous 
Red Galaxy (LRG) and Main samples, and also include the 2-degree Field Galaxy Redshift 
Survey (2dFGRS) data. In total, this sample comprises 893 319 galaxies over 9100deg 2 . 
Baryon Acoustic Oscillations are observed in power spectra measured for different slices 
in redshift; this allows us to constrain the distance-redshift relation at multiple epochs. We 
achieve a distance measure at redshift z = 0.275, of r s (z d )/D v (0.275) = 0.1390 ± 0.0037 
(2.7% accuracy), where r s (zd) is the comoving sound horizon at the baryon drag epoch, 
Dv(z) = [(1 + z) 2 D\cz/H(z)] 1 / 3 , Da(z) is the angular diameter distance and H{z) is 
the Hubble parameter. We find an almost independent constraint on the ratio of distances 
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Dy(0.35)/Dy(0.2) = 1.736 ± 0.065, which is consistent at the 1.1 a level with the best fit 
ACDM model obtained when combining our z — 0.275 distance constraint with the WMAP 
5-year data. The offset is similar to that found in previous analyses of the SDSS DR5 sample, 
but the discrepancy is now of lower significance, a change caused by a revised error anal- 
ysis and a change in the methodology adopted, as well as the addition of more data. Using 
WMAP5 constraints on Q^h 2 and ft c h 2 , and combining our BAO distance measurements with 
those from the Union Supernova sample, places a tight constraint on Cl m = 0.286 ± 0.018 and 
Hq = 68.2 ± 2.2 kms" 1 Mpc -1 that is robust to allowing ^ and w ^ — 1. This result is 
independent of the behaviour of dark energy at redshifts greater than those probed by the BAO 
and supernova measurements. Combining these data sets with the full WMAP5 likelihood 
constraints provides tight constraints on both = —0.006 ± 0.008 and w = —0.97 ± 0.10 
for a constant dark energy equation of state. 

Key words: cosmology: observations, distance scale, large-scale structure of Universe 



1 INTRODUCTION 

"What is the nature of dark energy?" is one of the current key ques- 
tions in physical science. Distinguishing between competing the- 
ories will only be achieved with precise measurements of the cos- 
mic expansion history and the growth of structure within it. Among 
current measurement techniques for the cosmic expansion, Baryon 
Acoustic Oscillatio ns (BAO) appear to ha ve the lowest level of sys- 
tematic uncertainty dAlbrecht et alj2006h . 

BAO are a series of peaks and troughs, with a wavelength of 
approximately 0.06/iMpc -1 that are present in the power spec- 
trum of matter fluctuations after the epoch of recombination, and on 
large-scales. They occur because the primordial cosmological per- 
turbations excitesOTnd_wavesjnthe relativ i stic plasma of the early 
universe dsilkl 1 19681; IPeebles & Yul |l97Ct ISunvaev & Zel'dovichl 
ll97d : lBond & Efstathioj|l984ll98llHoltzmanlll989h . Radiation 
pressure drives baryonic material away from the seed perturbations 
until the ionised material recombines at redshift z ~ 1000. The 
momentum of the baryonic material means that the motion con- 
tinues for a short time after recombination, until an epoch known 
as the baryon-drag epoch. The wavelength of the BAO is related 
to the comoving sound horizon at the baryon-drag epoch, which 
depends on the physical densities of matter Vl m h 2 and of baryons 



dKomatsu et alj2009T) give that r s (zh) — 153.5 Mpc (see Section[7] 
for details). 

BAO occur on relatively large scales, which are still pre- 
dominantly in the linear regime at present day; it is therefore 
expe cted that BAO should also be seen in the galaxy distribu- 



tion {Goldberg & Straus Jl9 98: Meiks in et aLHl999l;ISpringel et al.l 



l2005l ; ISeo & Eisensteirj|2005l ; IWhitejl2005l ;lEisenstei n et alj|2007l) . 
We can therefore use BAO as standard rulers to constrain the expan- 
sion of the Universe if the comoving sound horizon at the baryon 
drag epoch is known. The apparent size of the BAO measured 
from observations then leads to measure ments of the Hubble pa- 
rameter and the angular diameter distan ce ( Seo & Eisenstein 2003 ; 
iBond & Glazebrookil2003l ; lHu & HaimaJl2006r . lMatsubadl2004h . 

The aco ustic signature has now been convincingly detected at 
low r edshift (Per cival et al.l200"ll ; ICole et alj 2005; Eis enstein et al.l 
120051 ; lHuetsill2006h usin g the 2dF Galaxy Redshift Survey (2dF- 



GRS;IColless et alj2003l) and the Sloan Digital Sky Survey (SDSS; 
lYork et alj 12000). The detection has subsequently been refined 
using more data and better techniques, and is now producing 
compe titive constraints on cosmological models. iTegmark et al.l 
( 2006) analysed the Sloan Digital Sky Survey (SDSS) Data Release 



4 (DR4; Adelm ai>McCarthy et"al]|2006t) Luminous Red Galaxy 
(LRG) sample. iPercival et a l. (2007a. b) presented the power spec- 
trum of the Sloan Digital Sky Su rvey (SDSS) Data Release 5 (DR5; 
lAdelman-McCarthv et al.ll2007l) galaxy sample and considered the 
shape of the power sp ectrum and measured the matter density us- 
ing the BAO features. IPercival et alj d2007d ) took this analysis a 
stage further by fitting the SDSS data, combined with the 2dFGRS, 
with mod els of the dis t ance- redshift relation. iGaztanaga et al.1 
( 120081) andlSanchez et alj d2009l) have also analysed the SDSS DR6 
dAdelman-McCarthv et alj 12008) sample, obtaining cosmological 
constraints from th e radial and spheri cally averaged BAO signal. 
In a recent analysis. iKazin et aT] ( l2009h have calculated the correla- 
tion function of the SDSS DR7 dAbazaiian et"ai1 f2009) LRG sam- 
ple, and have shown that their results agree with those presented in 
our paper. Two studies have also considered the clustering of the 
LRGs at high redshift within the SDSS survey, using photometric 
redshifts to estim ate galaxy distances IPadmanabhan et al]|2007l ; 
iBlake et alj2007l) . 

In this paper, we analyse the clustering of galaxies in the 
spectroscopic SDSS DR7 sample, including both LRG and Main 
galaxy samples, combined with the 2dFGRS, and measure the 
BAO signal in a series of redshift slices. SDSS DR7 marks the 
final release of galaxies observed using the standard SDSS tar- 
geting algorithm, and the sample we analyse covers a solid an- 
gle of 7930 deg 2 , including a 7190 deg 2 contiguous region in the 
North Galac tic Cap. The Baryon Oscillation Spectroscopic Sur- 
vey (BOSS; ISchlegel et al.|[2009ah . part of the SDSS-III project, 
will adopt a different targeting algorithm, focusing on galaxies and 
quasars at higher redshifts. 

The observed amplitude of the large-scale galaxy cluster- 
ing d e pends on both galaxy colour and luminosity dTegmark et al.1 
120041 ; Izehavi et alj|2005l; ISwanson et alj|2008h . Using the SDSS 
DR5 sample, ICresswell & Percivall ([2009) showed that for blue 
galaxies, the deviation in the shape of the galaxy power spectrum 
from the linear matter power spectrum at k > 0.1 /iMpc -1 is a 
strong function of luminosity, while it is almost constant for red 
galaxies. It is therefore difficult to extract the underlying matter 
power spectrum from a galaxy power spectrum measured for a pop- 
ulation of galaxies where the distribution of galaxy colours and lu- 
minosities changes with spatial location, such as that provided by a 
magnitude-limited catalogue. In contrast, the luminous red galaxy 
population, which comprises the high redshift part of the sample 
analysed here, has a simpler relation with the m atter field, in tha t 
there is a single galaxy population to consider dReid et al]|2008h . 
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In a companion paper dReid et al]|2009h . we apply a grouping al- 
gorithm to recover the halo power spectrum from the LRGs, then 
calibrate the relation of the halo power spectrum to the linear the- 
ory power spectrum using simulations. We are then able to extract 
cosmological information from the large scale shape of the power 
spectrum in addition to the BAO signal, though the constraints are 
more tightly embedded in the assumed cosmological framework. 

BAO in the galaxy power spectrum are only weakly af- 
fected by the effects of non-linear structure formation and scale- 
dependent galaxy bias, because they are on such large scales. 
The primary consequence is a damping on small sc ales, which 
can be well approximated by a Gaussian smoothing dBharadwail 
19961: ICrocce & Scoccimarroll2006t [200^: lEisenstein et al.l I2067F 



Matsubara 2008a, b). The observed BAO, defined as the ratio of the 
observed power spectrum P b s to a smooth fit to this power P n w, 
BAOohs = -Pobs/^nw, are related to the original BAO in the linear 
matter power spectrum BAOn n , defined similarly, by 



BAOohs = Gdamp-BAOlin + (1 — Gda 



(1) 



where Gdamp = exp(— D damp ), and the damping scale, 
D damp is se t to 10 h~ Mpc for redshift-space power spectra 
at 2 ~ 0.3 teisenstein et al.ll2007t) . This damping of the linear 
power is a relatively benign effect as it does not affect the po- 
sitions of the BAO, although it does reduce the signal available. 
Additional, more pernicious effects such as the mixing of modes 
in the power spectrum, can generate shifts in the BAO position 
dCrocce & Scocci marro 200 8]); for biased tracers, these offsets can 
be at the percent level dSmifh et alj|2007t) . and are therefore impor- 
tant as we wish to make percent level distance measurements. 

In our analysis, we measure BAO relative to a model that al- 
lows for smooth changes in the underlying shape of the power spec- 
trum, which alleviates some of this shift. Physical models of BAO 
positions in observed re dshift-space power spectra rela t ive to a such 
a fitted smooth model dCrocce & Scoccimarrdl2008l; Smith et al. 
l2008l : ISanchez et al.ll2009l: IPadmanabhan & Whitell2009h, and nu- 
meric al simulations dAngulo et al.l2008l : ISeo et al.l2008l : lKim et ail 
2009) suggest we should expect residual shifts at the sub-percent 
level. These are below the precision of current experiments: e.g. 
in this paper we present a BAO distance scale measurement with 
2.7% accuracy. Therefore, we adopt a procedure that allows for 
the damping as well as smooth changes in the underlying shape of 
the power spectrum, but no more. The analysis of future surveys, 
which will lead to tighter distance-redshift constraints, will clearly 
also have to allow for non-linear effects, either by physical model- 
ing, simulations, or by usin g methods which attempt to reconstruct 
the initial fluctuation fiel d JEisenstein et alj|2007l : beo et alj|2008l ; 
IPadmanabhan et alll2009l) . 

The SDSS and 2dFGRS data are discussed in Sections|2]and[3] 
The basic methodolog y, presented in Section[4] is similar to that of 
IPercival et al.l (12007c). although we have revised the calculation of 
the window function to increase the computational speed. We also 
perform an extensive test of the derived errors, running mock cata- 
logues through our full analysis pipeline to test the confidence in- 
tervals quoted (Section|5j. Results are presented in Section[6]and[7] 
tested for robustness in Section[8]and placed in a cosmological con- 
text in Sections and [9] A comparison with our DR5 analyses is 
given in Section[l0]and we finish with a discussion in SectionfTTI 

In this paper we use the standard cosmological parameters. For 
flat ACDM models these are the Hubble constant Ho, the densities 
of baryonic matter ttb, cold dark matter tt c , all matter fl m , and dark 
energy Qa- Going beyond this simple class of models, we use the 
equation of state of the dark energy w, the curvature energy density 



Ok and total energy density f2tot- When combining with informa- 
tion from the CMB, we also consider some parameters that are not 
constrained by the BAO: r is the optical depth to re-ionization, n s 
is the scalar spectral index, and Aos is the amplitude of curvature 
perturbations at k — 0.05 Mpc . 



2 THE DATA 

The SDSS-I an d SDSS-II projects used a 2.5m telescope 
dGunn et al. l2006h. to obtain imaging data in fiv e passbands u, g, r, 
i and z dFukugita et alj| 19961: iGunn et alj|l998h. T he images were 
reduced dLupton et al.l200ll : IStoughton et al.l2002l;|Pier et alj2003l: 
Ivezic et al.N2004l) and calibrated dLupton et al .11 19991 ; iHogg et al.l 



200ll ; ISmith et al.ll2002l ; lTucker et alj|2006l) . and galaxies were se- 



lected i n two ways for foll ow-up spectroscopy. The main galaxy 
sample dStraussetai1l2002h targeted galaxies brighter than r 
1 7.77' (approximately 90 per square degree , with a weighted me- 
dian redshift z = 0.10). The DR7 sample dAbazaiian etai]|2009h 
used in our analysis includes 669 905 main galaxies dStrauss et al.l 
|2002|) with a median redshift of z = 0.12, selected to a limiting 
Galactic extinction-corrected Petrosian magnitude r < 17.77, or 
r < 17.5 in a small subset of the early data from the survey. The ef- 
fect of the inclusion of the early SDSS data is tested in Section [&2l 
In addit ion, our sample includ es 80 046 Luminous Red Galaxies 
(LRGs; lEisenstein et al.l 1200 lh . which form an extension of the 
SDSS spectroscopic survey to higher redshifts 0.2 < 2 < 0.5. 
Of the main galaxies, 30 530 are also classified as LRGs and are 
intrinsically luminous with Mo.i r < —21.8, where Mo.i r is the 
Galactic extinction and K-corrected r-band absolute galaxy magni- 
tude. We apply this requirement to all of our LRGs, so our sample 
includes 110 576 LRGs in total, with a weighted median redshift 
of z = 0.31. Although the main galaxy sample contains signifi- 
cantly more galaxies than the LRG sample, the LRG sample covers 
more volum e. Redshift distribution s for these samples are shown 
in figure 2 of lPercival etail d2007bh . In our default analysis we use 
SDSS P etrosian magnitudes calibr ated using the "uber-calibration" 
method dPadmana bhan et al 2008), although we test against data 
calcu lated using the original calibration methodology dTucker et al.l 
2006). Where specified, we have K-corre cted the galaxy luminosi- 
ties using the methodology outlined by Blant onet al.l J2003al lbh. 
Further details of th e cuts applied to the data can be found in 
|Percivaletai1d2007bl) . 

Due to the finite size of the fibers, spectra cannot be obtained 
for both objects in a p air closer than 55 arcs ec, within a single spec- 
troscopic tile. Tiling dBlanton et al.ll2003ah deals with this to some 
extent by allowing plate overlaps to provide multiple observations 
of crowded regions. Even so, not all galaxies in such r egions which 
meet t he target selection criteria could be observed. IZehavi et al.l 
d2002l) corrected for this undersampling by assigning the redshift 
of the nearest observed galaxy to a galaxy which was not observed 
due to crowding, and showed that this provides sufficient correc- 
tion for large-scale structure studies. We apply this correction in 
the present work, and test it to show that our results are insensitive 
to this in Section [8~Tl 

In order to increase the volume covered at redshift 2 < 0.3, we 
include 143 368 galaxies from the 2dFGRS sample. These galax- 
ies, selected to an ext inction-corrected m agnitude limit of approxi- 
mately bj = 19.45 (Colless et al. 20031) from regions of sky not 
covered by the SDSS sample, cover two contiguous regions to- 
talling ~ 1 200deg 2 . They do not include the 2dFGRS random fields, 
a set of 99 random 2 degree fields spread over the full south- 
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0.5 


129045 


0.27 


1.92 


7 


0.3 


0.5 


68 074 


0.15 


0.67 



V cS = / d\ 



n{r)P 
1 + n(r)P 



(2) 



Table 1. Parameters of the redshift intervals analysed. V e s is given in 
units of h~ 3 Gpc 3 , and was calculated as in Eq. j2) using an effective 
power spectrum amplitude of P = 10 4 /i — 3 Mpc, appropriate on scales 
k ~ 0.15 h Mpc -1 for a population with bias b = 1 . 7. The average galaxy 
number density in each bin n is in units of 10~ 4 ( h" 1 Mpc) 3 . 



ern Galactic cap, as these would complicate the window function. 
The galaxies cover < z < 0.3, with a weighted median at 
z = 0.17. The redshi ft distribution of the sample was analysed 
as in ICole et al.l d2005f) for < z < 0.3, and we use the same 
synthetic catalogues to model the unclustered expected galaxy dis- 
tribution within the reduced sample. 

We assume that each galaxy is biased with a linear determin- 
istic bias model, and that t his bias depends on Mo i r according to 
iTegmark et"ai] d2004h and IZehavi et all d2005h . All galaxies were 
weighted using this model so the fluctuation amplitudes match 
those of L, galaxies, where L* was calculated separately for the 
SDSS and 2dFGRS. We include an extra normalisation factor to 
the 2dFGRS galaxy bias model to correct the relative bias of L* 
galaxies in the different surveys. This was calculated by matching 
the normalisation of the 2dFGRS and SDSS bias-corrected power 
spectra for k < 0.1 ftMpc -1 . In principle, we could have added 
information on galaxy bias from the BAO, since the small-scale 
damping (see Eq. QJ depends on how strongly nonlinear the un- 
derlying dark matter density fluctuations are. As we show in Sec- 
tion [8T6J this information is limited for the current data, but future 
surveys may be able to exploit changes in this damping as a func- 
tion of galaxy properties, such as colour and luminosity. 



3 SPLITTING INTO SUB-SAMPLES 

In order to probe the distance-redshift relation in detail, ideally we 
would analyze BAO measured in many independent redshift slices. 
However, if the slices are too narrow in redshift, then there is in- 
sufficient signal and the BAO cannot be recovered with sufficient 
accuracy to give a likelihood with close to a Gaussian distribution 
(see the discussion in Section|5j- If the slices are too wide, or too 
many overlapping slices are chosen, the covariance matrix becomes 
close to singular, potentially leading to numerical instability. In or- 
der to balance these competing requirements, we have chosen to 
analyse the redshift slices presented in Table Q] The power spec- 
tra will be correlated, and these correlations, together with corre- 
lations of P(k) values at different k within each redshift slice, are 
included in the covariance matrices in our analysis. Note that we 
include slice 7, for which the effective volume is relatively small, 
because of the interesting redshift range covered. 

As well as giving the redshift limits of the slices in Table Q] 
we also give the number of galaxies in each including both the 2dF- 
GRS an d the SDSS, and the effective volume, calculated from the 
integral dFeldmanetaill 19941) 



where ra(r) is the observed comoving number density of the sam- 
ple at location r and P is the expected power spectrum amplitude. 
To calculate V a g for our redshift slices, distances were calculated 
assuming a fiducial flat ACDM cosmology with fl m = 0.25. For 
the numbers given in Table Q] we fix P = 10 4 /t" 3 Mpc 3 , appro- 
priate on scales k ~ 0. 15 h Mpc" 1 for a pop ulation with bias 
b = 1.7. For comparison. lEisenstein et ail d2005) analyse a sample 
with VcS = 0.13 /i _3 Gpc 3 , approximately a third of the effective 
volume of slice 1 . 

We fit models to three sets of power spectra: 

(i) We fit a single power spectrum for the SDSS LRG sample 
covering 0.15 < z < 0.5. 

(ii) We fit three power spectra for slices 1 , 3 and 6 approximatel y 
corresponding to the procedure adopted bv lPercival et alj d2007ch . 
Although we now use slices constrained by redshift rather than 
galaxy type, the < z < 0.3 slice is dominated by SDSS main 
galaxies, while the 0.2 < z < 0.5 slice is dominated by LRGs. 

(iii) We fit six power spectra for slices 2— >7, which allows a test 
of the distance-redshift relation at greater resolution. 

We consid e r opt ion [(1)1 to tie in with the analysis presented by 
iReid et all d2009l) . and to demonstrate the effect of collapsing 
the clusters in redshift-space where we try to reconstruct the 
halo power spectrum . Option |(ii)| is close to the approach of 
IPercival et aT] d2007ch . where the SDSS main galaxies and 2dFGRS 
galaxies were analysed separately from the SDSS LRGs. Option 
|(iii)| allows us to see if there is more information available beyond 
measuring the distance-redshift relation at two redshifts. The slices 
do overlap in redshift, but we will properly take into account the 
covariance between the results when we fit to cosmological param- 
eters. 



4 BASIC METHODOLOGY 

Power spect ra were calcula t ed for each catalogue usin g the Fourier 
method of iFeldman et al] dl994) . as applied by IPercival et al.l 
(2007b). In this method a weighted galaxy over-density field is de- 
fined and Fourier transformed, then the spherically averaged power 
is measure d. We use the l umino sity dependent galaxy weights ad- 
vocated bv lPercival et alj J2004I) . as described in Section|2] To con- 
struct the over-density field, we need to quantify the expected spa- 
tial distribution of galaxies, in the absence of clustering. The stan- 
dard method for this is to use an unclustered random catalogue, 
which matches the galaxy selection. To calculate this random cata- 
logue, we fitted the redshift d istributions of the galaxy samples with 
a spline fit dPress et al.ll 992"). and the an gular mask was dete rmined 
using a routine based on a H E ALPIX dGorski et al.ll2005l) equal- 
area pix elization of the sphere dPercival et alj2007bh . IPercival et al.l 
(2007b) used a random catalogue containing ten times as many 
points as galaxies. For the sparse LRGs, this approach induces 
significant shot noise, so we now use one hundred times as many 
random points as LRGs. We have also increased the resolution at 
which the radial distributio n of galaxies is quantified, now using a 
spline fit dPresset all 1992b with nodes separated by Az — 0.0025. 
As an alternative to this radial selection, we could have simply 
adopted the redshift of a randomly chosen galaxy for each of our 
points in the random catalogue. In Section [8?7l we show that these 
two possibilities give consistent results. 
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Galaxy redshifts were converted to distances using a fidu- 
cial cosmology (flat ACDM model with Q m — 0.25). For each 
distance-redshift model to be tested, we do not recalculate the 
power spectrum, but instead change the interpretation of the power 
spectrum computed assuming the fiducial ACDM galaxy distances. 
We do this through a window function, which relates the true and 
measured power spec tra. This follows the procedure adopted by 
iPercival et all d2007ch . but we now use a revised, computationally 
less intensive method for calculating the windows, as described in 
Appendix lAl 

A model of the BAO was created by fi tting a linear mat - 
ter power spectrum, calculated using CAMB dLewis et al.l 2000). 
which numerically solves the Boltzman equation describing the 
physical processes in the Universe before the baryon-drag epoch, 
with a cubic spline to remove the broad shape of the power, leav- 
ing the oscillations. The theoretical BAO were then damped with 
a Gaussian model as in Eq. (Q3, following the simulation results 
of lEisenstein et al] d2007l) . For our def ault fits, we assume tha t the 
damping scale Ddamp = 10 ft -1 Mpc dEisenstein et alj|2007l) . but 
we also consider fits where this scale is varied (Section |8.6| l. As 
discussed in Section Q] we do not attempt to correct for any shift 
induced by non-linear physics, because they are expected to be at a 
level below our statistical error. 

The power spectrum measured from the data was fitted by a 
model constructed by m ultiplying this BAO model with a cubic 
spline dPress et al J 1992b . which enables the power spectrum model 
to match the overall shape of the data power spectrum. Each power 
spectrum model was then convolved with a window function that 
corrects for both the survey geometry and the differences between 
our fiducial cosmological model used to convert redshift to dis- 
tances and the cosmological model to be tested (see Appendix lAl. 
The free parameters of the model are the nine nodes of the cubic 
spline fixed empirically at k = 0.001, and 0.025 < k «S 0.375 
with Afc = 0.05, and the parametrisation of Dy{z) used to cal- 
culate the correct window function. The spline nodes were refitted 
for every cosmology (or Dv(z)) tested. A power spectrum model 
with this spline n ode separation was te sted by fitting many mock 
power spectra by IPercival et alj d2007al) and was shown to match 
these without leaving significant residuals in the measured "shift" 
between BAO in the m odel and data power spectra. This approach 
was also considered bv lSanchez et"al| {2008). who found that it did 
not induce a bias in the recovered BAO constraints. 

For a redshift survey in a thin shell, the position of the BAO 
approximately constrains d z = r s (zd)/Dy(z), where r a (zd) is 
the comoving sound horizon at the baryon drag ep o ch, Dy (z) = 
[(1 + z) 2 D 2 A cz/H(z)] 1/3 dEisenstein etai]|2005l : IPercival et all 
l2007ch . Da is the angular diameter distance, and H(z) is the Hub- 
ble parameter. We see that, although our power spectrum fitting 
procedure measures Dy (z) for a fixed BAO model, we should con- 
sider the constraints as measurements of d z , with r s (zd) calculated 
for the flat ACDM model for which we created the BAO model, 
r s (zd) = 111.4 h~ M pc = 154.7 Mpc, using equation 6 of 
lEisenstein &Hu]d 19981) . and assuming h = 0.72, Q. b h 2 = 0.0223, 
and Q m — 0.25. This value of r s (zd) is only used to index this 
model: as described above, the actual BAO model was calculated 
from a power spectrum predicted by CAMB. If the constraints pro- 
vided in this paper are to be used to constrain a set of models where 
r a (zd) for this fiducia l model is calculated in a different way (i.e. 
not using equation 6 of Eis enstein & Hul 19981) , then our constraints 
should be adjusted to match. 

The comoving distance-redshift relation is modelled as a cu- 
bic spline in the parameter Dy(z). We consider models for Dy(z) 



k / h Mpc 1 
0.02 0.05 0.1 0.2 




-1.5 -1 -0.5 

log 10 k / h _1 Mpc 



Figure 1. Average power spectra recovered from the Log-Normal cata- 
logues (solid lines) compared with the data power spectra (solid circles with 
1-<t errors) for the six samples in Table[TJ The errors on the data were calcu- 
lated from the diagonal elements of the covariance matrix calculated from 
these log-normal catalogues. The power spectra have been offset by 0.5dex 
for clarity, with the upper power spectrum having the correct normalisation. 



with two nodes at z — 0.2 and z = 0.35, or with four nodes at 
z = 0.1, z = 0.2, z = 0.3, z = 0.45. Results are presented as 
constraints on d z . The error between cubic spline fits to Dv(z) 
with two nodes at z = 0.2 and z = 0.35, t o the ACDM distance - 
redshift relations was shown in figure 1 of IPercival etafl d2007cl) . 
and is < 1% for a flat ACDM cosmology with fi m = 0.25 at 
0.15. 

Power spectra are presented for the redshift slices described in 
Section[3]in Fig.Q] for 70 band powers equally spaced in 0.02 < 
k < 0.3/i Mpc" 1 . We see that the power spectra from the different 
redshift intervals are remarkably consistent, with P(k) decreasing 
almost monotonically to small scales. 

In order to calculate the covariances between th e data, we have 
created 10000 Log-N ormal (LN) density fields dColes & Jonesl 
ll99ltlCole et alj2005h from which we have drawn overlapping cat- 
alogues for each of our 7 redshift slices. Catalogues were calculated 
on a (512) 3 grid with box length 4000 h' 1 Mpc. Unlike N-body 
simulations, these mock catalogues do not model the growth of 
structure, but instead return a density field with a log-normal distri- 
bution, similar to that seen in the real data. The window functions 
for these catalogues were matched to that of the 2dFGRS+SDSS 
catalogue with the original calibration. The input power spectrum 
was a cubic spline fit matched to the data power spectra (i.e. the 
smooth part of our standard model), multiplied by our default 
damp ed ACDM BAO model calculated using CAMB dLewis et all 
2000). The LN power spectra were used to determine a covariance 
matrix between slices and for different band powers in each slice, 
assuming that the band-powers were drawn from a multi-variate 
Gaussian distribution. Average recovered power spectra for each 
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redshift interval are compared with the data power spectra in Fig.Q] 
Clearly the general shape of the average power spectra of the LN 
catalogues is well matched to that recovered from the data. Using 
the inverse of this covariance matrix, we estimate the likelihood 
of each model assuming that the power spectra band-powers for 
0.02 < k < 0.3 /iMpc -1 were drawn from a multi-variate Gaus- 
sian distribution. 



5 TESTING THE ANALYSIS METHOD WITH MOCK 
DATA 

5.1 The model fit 

We now consider using a subset of our LN catalogues to test 
our analysis procedure. For 1000 of the mock catalogues, we fit 
spline xBAO models to extract distance constraints from the BAO, 
as described in Section|4] A small average shift of 1.3% in the BAO 
scale was recovered between the power recovered from the LN cat- 
alogues, and the input power spectrum used to create them. If we 
correct the 1000 power spectra measured from the LN catalogues 
by multiplying each power spectrum by the expected power divided 
by the average recovered power spectrum, the average shift drops 
below 0.3%, well within la. 

To test the origin of the observed 1.3% shift, we have also 
drawn 1000 power spectrum realisations from a multi-variate Gaus- 
sian distribution with covariance and mean matched to those of the 
data. These mock catalogues were fitted using the procedure de- 
scribed in Section[4] No shift in the BAO position was found from 
the fits to these catalogues, within the statistical limits of the anal- 
ysis (~ 0.3%). The distribution of recovered distance constraints 
was well matched to that recovered from fitting the corrected LN 
power spectra. Thus the 1.3% shift described above must be due 
to the LN procedure itself. The expected shift is dependent on the 
statistic used to measure the BAO position. The Log-Normal corre- 
lation function £,ln, and Gaussian correlation function £a{r) of a 
field with the same power spectrum but with Gaussian statistics, are 
related by 1 + £ljv(V) = exp[^c{r)]. If we na d used the peak in 
the correlation function as our standard ruler then, for the LN cata- 
logues, we would have expected no BAO shift. However, the same 
is not true of our BAO x spline model fitting procedure, which fits 
the BAO in the power spectrum over a range of scales. 

Numerical simulations offer a better way to model the true 
Universe, and recent results from simulations show that we should 
expect a less significant shift between the BAO positions in the lin- 
ear matter and gala xy power spectra than th e 1.3% shift found for 



the LN catalogues ( Seo & Eisenstein 



ISeo & Eise nstein 2007; Anguloetal 



ISpringel et alj|2005l : 

2008). The exact shift re- 



2003 



quired for the catalogues we analyse is not well constrained by 
these simulation results, and we consequently do not alter our anal- 
ysis to include such a shift. 



5.2 The likelihood surface 

We use the Gaussian and LN power spectra samples to assess the 
nature of the likelihood for the BAO scale recovery. We consider 
fits to either three or six power spectra as described in Section fj] 
parametrising Dv{z) with a cubic spline with two non-zero nodes 
at z — 0.2 and z — 0.35. For each of the 1000 fits, we have mea- 
sured the difference between the maximum likelihood value and 
the likelihood at the parameters of the true cosmological model. 



-21n£/£ truc 




fraction of samples 






3 slices 


6 slices 




standard 


revised 


standard 


revised 


<2.3 


0.579 


0.666 


0.551 


0.667 


<6.() 


0.892 


0.946 


0.862 


0.948 


<9.3 


0.966 


0.983 


0.955 


0.981 



Table 2. Fraction of fits to the Log-Normal power spectra in which the ratio 
of the likelihood maximum and the likelihood for the true cosmological 
model is less than the given limit. For a Gaussian likelihood, these limits 
correspond to 68%, 95% and 99% confidence intervals. We show results 
where we have corrected the errors as described in the text by multiplying 
the band-power errors by 1.14 for three redshift slices, and 1.21 for six 
redshift slices. 



The fraction of samples with — 21n£/£ t ruc < 2.3, 6.0, 9.3, cor- 
responding to 68%, 95% and 99% confidence intervals, are given 
in Table [2] We find that in order to match the expected numbers 
of samples with likelihoods within the standard lcr Gaussian confi- 
dence intervals, we must increase the errors on the power spectrum 
band powers by 14 ± 2% if we fit to three power spectra. For fits to 
six power spectra, we must increase the errors by 21 ± 2% in order 
to match the expected la Gaussian confidence intervals. Although 
in this paper we do not consider fitting to a single power spectrum, 
we have repeated this analysis for BAO fits to the LRG sample of 
iReid et all (2009), and find that we must increase the errors on the 
power spectrum band powers by 10 ± 2% to match the expected 
confidence intervals. 

Because the same increase in the confidence intervals is re- 
quired for both LN and Gaussian mock catalogues, this change 
must be caused by the methodology of fitting BAO, rather than the 
Gaussian to Log-Normal density field transition. In fact, we believe 
that it is caused by the non-Gaussian nature of the likelihood sur- 
face. We should expect the likelihood surface to be non-Gaussian 
to some extent in any case because there is a minimum in the likeli- 
hood where the observed and model BAO are perfectly out of phase 
in A;-space: this represents the worst possible match between data 
and model. Adjusting the covariance matrix to match the distri- 
bution of best-fit distance-scales to the expected 68% confidence 
interval does not quite match the 95% or 99% confidence intervals, 
although it corrects for most of the difference. This shows that the 
confidence intervals cannot perfectly match those for a Gaussian 
distribution. 

To test this further, we have created a set of 1000 Gaussian 
power spectrum realisations with errors that are 10% of those in our 
standard sample. For these catalogues, the distribution of best-fit 
Dy{z) matches that expected from the likelihood distribution un- 
der Gaussian assumptions. No correction is required, and the likeli- 
hood distribution is much closer to that for a multi-variate Gaussian 
distribution around the likelihood maximum. Thus the requirement 
to increase the errors on the data disappears when we fit less noisy 
data, as we would expect if it is caused by fitting noisy data, which 
is giving a non-Gaussian likelihood surface. 

The average likelihood surfaces measured from our 1000 fits 
to sets of three power spectra and six power spectra drawn from LN 
catalogues are shown in Fig. [2] We also plot the centre of the local 
likelihood maxima nearest to the input cosmological parameters for 
each model. The fractions of points within each contour are given 
in Table [2] the errors on the power spectrum band powers have 
been adjusted for each plot as described above so that ~68% of the 
points lie within the —2 In £ = 2.3 contour. 
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Figure 2. Average likelihood contours recovered from the analysis of three 
power spectra (top panel) and six power spectra (bottom panel) measured 
from 1000 Log-Normal density fields. Contours are plotted for —2 In C = 
2.3, 6.0, 9.2, corresponding to two-parameter confidence of 68%, 95% and 
99% for a Gaussian distribution. Contours were calculated after increas- 
ing the errors on the power spectrum band-powers as described in the text. 
Solid circles mark the locations of the likelihood maxima closest to the 
true cosmology. We have plotted the likelihood surface as a function of 
Dy(2)/Mpc, for fixed r s (z,i) = 154.7 Mpc, to show distance errors if 
the comoving sound horizon is known perfectly. The values of Dy for our 
input cosmology are shown by the vertical and horizontal solid lines. 



6 RESULTS 

Baryon Acoustic Oscillations are observed in the power spectra 
recovered from all redshift slices of the SDSS+2dFGRS sample 
described in Section [3] and are shown in Fig. [3] where we plot 
the measured power spectra divided by the spline component of 
the best-fit model. In our default analysis we fit power spectra 
from six redshift slices as described in Section [3] using a spline 
for Dy(z) with two nodes at z — 0.2 and z = 0.35. We as- 



sume a fixed BAO damping scale of Dda 
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fit to all SDSS and non-overlapping 2dFGRS data. The effect of 



0.1 0.2 
k / h Mpc -1 

Figure 3. BAO recovered from the data for each of the redshifts slices 
(solid circles with 1-cr errors). These are compared with BAO in our de- 
fault ACDM model (solid lines). 



these assumptions is considered in Section [8] The resulting like- 
lihood surface is shown in Fig. [4] as a function of Dy (z) /Mpc, 
for fixed r 3 (zd) = 154.7 Mpc, to show distance errors if the co- 
moving sound horizon is known perfectly. The constraints should 
be considered measurements of r s (zd)/Dv(z) (see Section |4j. 
Fig. [4] reveals a dominant likelihood maximum close to the pa- 
rameters of a ACDM cosmology with Q m = 0.25, h = 0.72, & 



© 0000 RAS, MNRAS 000, 000-000 



8 W.J. Percival et al. 



o 
a 
2 



LO 

CO 



Q 



o 
o 



o 
o 
o 




600 



800 



1000 



D v (z=0.2) / Mpc 



Figure 4. Likelihood contour plots for fits of two Dy(z) cubic spline 
nodes at z = 0.2 and z = 0.35, calculated for our default analysis 
using six power spectra, uber-calibration, a fixed BAO damping scale of 
D damp = 10 h,- 1 Mpc, and for all SDSS and non-overlapping 2dFGRS 
data. Solid contours are plotted for —2 In £/£truc < 2.3, 6.0, 9.3, which 
for a multi-variate Gaussian distribution with two degrees of freedom cor- 
respond to 68%, 95% and 99% confidence intervals. Likelihoods were ad- 
justed to match these Gaussian confidence intervals as described in Sec- 
tion|3] We have plotted the likelihood surface as a function of Dy (#) /Mpc, 
for fixed r.« (z,/) = 154.7 Mpc, to show distance errors if the comoving 
sound horizon is known perfectly. We also show a multi-variate Gaussian 
fit to this likelihood surface (dashed contours). The values of Dy for a flat 
ACDM cosmology with Q m = 0.25, h = 0.72, & Q, b h 2 = 0.0223 are 
shown by the vertical and horizontal solid lines. 



Q b h 2 = 0.0223. There are also weaker secondary maxima at lower 
£V(0.2), which are considered further in Section [878l The signif- 
icance of detection of BAO corresponds to Ax = 13.1, which 
is approximately 3.6a. As this is relative to an arbitrary smooth 
model, this test is more general, and hence the s ignificance can- 
not be directly compared with results presented by Eisenstein et al. 
fe005h . 

We have matched the likelihood surface shown in Fig. [4] 
around the dominant maximum to a multi-variate Gaussian model. 
Using this Gaussian fit, we find that the best fit model has 



d .2 = 0.1905 ±0.0061 (3.2%), 
do.35 = 0.1097 ±0.0036 (3.3%), 



(3) 



where d z = r s (zd) /Dy(z). These results are correlated with 
correlation coefficient r = 0.337. For a cosmological distance- 
redshift model with d z the likelihood can be well approximated by 
a multi-variate Gaussian with covariance matrix 



C : 



(Ado. 2 Ado. 2 ) 

(Ad .35Ad .2) 



(Ado.2Ado.35} 
(Ado.35Ado.35) , 



where Ad z = d z — d z . C has inverse 
1 / 30124 -17227 



C 



-17227 86977 



(4) 



(5) 



Without correcting the covariance matrix using the results 
from fitting to the LN power spectra as described in Section[5] the 
original average errors on do. 2, and do. 35 were 0.0051 and 0.0029, 
16% and 24% lower than those in Eq. (|3j, respectively. Compare 
with the band power errors which were increased by 21%, and we 
see that there is not a direct relation between changes in the band 
power errors and errors on d z , because of the non-linear nature of 
the fit. 

We diagonalise the covariance matrix of do. 2 and do. 35 to get 
quantities x and y 



1 



1.76 
1.67 



do. 2 

do. 35 



which gives 



x = 0.3836 ± 0.0102 
y = -0.0073 ± 0.0070. 

The distance ratio / = ZV(0.35)/.Dv(0.2) is given by 
1.67- \.76y/x 



f 



1 + y/x 



1.67 



■ My, 



(6) 



(7) 
(8) 



(9) 



where the last approximation neglects the small variations around 
the best-fit value of x = 0.3836; these would come to 0.002 in 
/, which is well within the errors. Thus, a; is a measurement of 
distance for the concordance cosmology and y is the deviation from 
the concordance distance ratio: x is measured to about 2.7%. y is 
consistent with zero to within about la. 

To high accuracy, the constraint x can be written as a con- 
straint on the distance to some redshift 0.2 < z < 0.35. In fact, 
r s (zd) / Dy (0.275) predicts x = do. 2 ± 1.76do.35 to a peak-to- 
peak precision of 0.04% over the range 0.05 < f2 ra < 1 (assuming 
a flat cosmology with w = —1). Thus, we can quote the x mea- 
surement as a measurement of do. 275 and quote the y measurement 
as a statistically independent measure of /. 

For the best-fit solution we have do. 275 = 0.362a;, giving 



do. 2 



0.1390 ±0.0037(2.7%). 



(10) 



We also have the statistically independent constraint 

/ = £V(0.35)/Dv(0.2) = 1.736 ± 0.065. (11) 

/ = 1.67 for our ACDM concordance cosmology, while SCDM 
with Q m — 1, J1a = has / = 1.55, which is only 2.9a from 
this result. Our constraint from the distance ratio only separates the 
concordance model from Q. m — 1 at 1.8<r, i.e., it is not a strong 
cosmological constraint, compared with the constraint on do. 275- 



7 COSMOLOGICAL INTERPRETATION 

We now consider how our constraints can be mapped into the 
standard basis of cosm ological parameters. From equation 6 of 
lEisenstein & Hul Jl998h . the sound h orizon can be approx imated, 
around the WMAP5 best-fit location dKomatsu et alj|2009h as 



r s (zd) 



153.5 



0.02273/ \ 0.1326 

Setting r 3 jid = 153.5 Mpc, and using Eq. (Hi we have 
ZM0.275) = (lW4±30)\r s (z d )/r sJid (z d )] Mpc 

^1 m h 



Mpc. (12) 



0.1326 



Mpc, (13) 
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and / = 1.736 ± 0.065 as our two statistically independent con- 
straints. 

The constraint on Dy (0.275), combined with a measurement 
of n m h 2 fromWMA P5 dffinshaw et alj2009l ; lDunklev et al.l2009l ; 
iKomatsu et aii r2009), is enough to measure fi m and Hq given in- 
formation about the distance scale from z = to z — 0.275. If the 
distance measure were at z = 0, then we would have a standard 
ruler defined by the CMB with which we could measure Ho, and 
combining this with f2 m /i 2 would yield Q m . In practice, one has to 
include the small corrections to Dy (0.275) that arise from the low- 
redshift cosmology. Noting that LV(0.275) = 757.4 h' 1 Mpc 
for a flat Q m — 0.282 ACDM cosmology, we can write h = 
\/{l m h 2 /\/fl m , and solve 



Q m = (0.282 ±0.015) 



0.1326 

Dv(z = 0.275, n m = 0.282) 
D v (z = 0.275) 



(14) 



where we have dropped the dependence of the sound horizon on 
Slbh 2 , which the WMAP5 data already constrains to 0.5%, 5 times 
below our statistical error. 

We can perturb the ratio of distances around the best-fit Q m = 
0.282, to give 

D v {z = 0.275) 



Dy(z = 0.275, On = 0.282) 



0.282 



[1 - 0.108fi* - 0.099(1 + w)] . (15) 



Using this approximation, we can manipulate Eq. d 1 4b to give con- 
straints on either Q m or h 



(0.282 ±0.018) 



0.1326 y 

x [1 ±0.25fi fc + 0.23(1 + iw)] . 

0.21 

h = (0.686 + 0.022) 



(16) 



0.1326 / 

X [1 - 0.13f} fe - 0.12(1 + w)] . (17) 

The additional uncertainty in Q m , ±0.018 in Eq. d!6b compared 
with ±0.15 in Eq. dl4b , is produced by the dependence of the dis- 
tance ratio on f2 m , In Eqns. d 1 6b & l !17b , the uncertainty in the first 
terms are correlated so as to leave Q m h 2 constant. One should ad- 
ditionally include the errors from Q, m h 2 , 0,^, and w, although these 
are consistent between the two results. 

Looking at the fractional error in Q m , the contribution from 
the uncertainty in the SDSS acoustic scale is about 6%, that from 
the uncertainty in fl m h 2 is about 2%, that from w is about 3% if 
the error on w is 10%, and that from curvature is below 1% unless 
the cosmology is rather non-standard. Hence our result is still lim- 
ited by the SDSS-II BAO data volume and not by our knowledge 
of the other cosmological parameters in Eq. dl6b . Of course, these 
expressions only hold for mild perturbations from the concordance 
cosmology; for other cases, one should return to the raw distance 
constraints. We note that these expressions have not used the angu- 
lar acoustic scale in the CMB, so they are independent of what is 
happening with dark energy at z > 0.35. 

Fig. [5] shows the BAO constraints from Eq. dl3b on fi m and 
SIa for ACDM cosmologies (upper panel), and on Q, m and w for 
flat models where constant w 7^ —1 is allowed (lower panel). 
We take a Gaussian prior of flmh 2 = 0.1326 ± 0.0063 and as- 
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Figure 5. Cosmological constraints on ACDM cosmologies (upper 
panel) and flat CDM models where we allow w to vary (lower 
panel), from WMAP5 (blue), Union supernova (green) and our con- 
straint on r s /£V(0.275) (solid contours). Contours are plotted for 
— 21n£/£truc < 2.3, 6.0, corresponding to 68% and 95% confidence 
intervals. The dashed lines show flat models (upper panel) and A models 
(lower panel). 



sume that the error on flb h 2 is negligible as th e WMAP5 data al- 
ready constrain it to 0.5% dKomatsu et al.ll200^) . These constraints 
exclude the angular acoustic scale in the CMB, so they are in- 
dependent of the dark energy behaviour at the redshifts beyond 
our sample. For com parison we plot the full WMAP5 constraints 
dKomatsu et alj2 009). which include the constraints on the distance 
to last scattering, and constraints from the Union supernova sample 
dKowalski e t al. 2008), which constrain angular diameter distance 
ratios up to z ~ 1. Results from full likelihood fits combining these 
data are presented in Section|9] 
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do. 2 cfo.35 

(a) default 0.1905 ± 0.0061 0.1097 ± 0.0036 

(b) no early SDSS, 2dFGRS 0.1923 ± 0.0072 0.1102 ± 0.0041 

(c) no2dFGRS 0.1907 ± 0.0062 0.1090 ± 0.0036 

(d) no early SDSS 0.1917 ± 0.0069 0.1109 ± 0.0044 

(e) fit to three P(k) 0.1901 ± 0.0066 0.1080 ± 0.0043 

(f) original calibration 0.1919 ± 0.0071 0.1094 ± 0.0046 

(g) varying D damp 0.1918 ± 0.0080 0.1100 ± 0.0048 

(h) (n(«)) sampling galaxies 0.1890 ± 0.0068 0.1102 ± 0.0045 

Table 3. Measurements of d z = r s (z^)/ Dy(z) at z = 0.2 and z = 0.35 
from the different analysis runs described in the captions to Figs.[7]&[H 



0.1 



0.2 



h Mpc 



Figure 6. The log ratio between the BAO recovered from the SDSS LRG 
power spectrum and the power spectru m of the halo ca talogue derived 
from the LRG sample as described by iReid et al (solid circles). 

For comparison we plot the BAO expected for a flat ACDM model with 
tt m = 0.25, h = 0.72, & U b h 2 = 0.0223 (solid line), and the errors on 
each measurement (grey shaded region). There are no oscillatory features 
induced by the cluster-collapse procedure, and the scatter is well within the 
errors. 



8 TESTING THE ROBUSTNESS OF THE RESULTS 

8.1 The effect of redshift-space distortions 

We have fitted our spline x BAO mo del to the observe d SDSS 
LRG power spectrum, as calculated by IReid et al ] d2009l) . where 
the galaxy power spectrum and derived cosmological constraints 
ar e presented . Using numerical simulations, a scheme is presented 
in lReid et al.l d2009 ) to recover the halo power spectrum from the 
LRG distribution by only keeping a single LRG within each halo. 
We have fitted both the galaxy and the halo power spectra with our 
spline x BAO model. The log ratio between the BAO recovered 
in the resulting fits is shown in Fig. [6] This shows that the cluster- 
collapse correction for these galaxies results in a smooth change 
in the power spectrum on the scales fitted, and does not alter the 
position or amplitude of the BAO in a significant way. 

Because of the different galaxy properties within the SDSS 
main galaxy sample, and the 2dFGRS, we do not attempt to correct 
for the more complicated distribution of galaxies within the haloes 
of that sample, and recover the halo power spectrum. In contrast, 
the halo occupation distribution of the SDSS LRGs is simple, in 
that there is only a single population of galaxies that are predom- 
inantly central ra ther than satellite galaxies in their hosting haloes 
dReid et alj|2009h . But we have seen that for LRGs, the correction 
is smooth, and we expect this to be true for the galaxies at z < 0.2 
as well. 



8.2 Sample selection 

We have run our full analysis pipeline using three subsamples of 
galaxies. Results from fits to Dy(z) with two nodes are shown in 
Fig- 13 for different catalogues, given r s (zd) = 154.7 Mpc. The 
best-fit constraints for these models on d z are given in Table|3] Our 
default analysis is included in panel (a) for comparison. Here, we 
analyse data from the SDSS and the 2dFGRS, including the early 
SDSS data, where we cut the sample at the extinction-corrected 



magnitude limit r < 17.5. We compare with results obtained (b) 
excluding the early SDSS data and the 2dFGRS, (c) using just the 
SDSS data, and (d) excluding the early SDSS data but including the 
2dFGRS. Including the early SDSS galaxies decreases the errors at 
redshift z = 0.2 and z = 0.35 by approximately 14%. Including 
the 2dFGRS galaxies has a smaller effect, decreasing the error at 
z = 0.2 by approximately 4%. The parameters of the best-fit so- 
lutions do not move significantly with any of the sample changes: 
do. 2 moves by a maximum of 0.3<r, while do. 35 moves by a max- 
imum of 0.2a. The inclusion of the 2dFGRS actually moves the 
best-fit solution for D v (0.35) / D v (0.2) slightly towards that of a 
concordance ACDM model. 



8.3 The number of redshift slices included 

We now consider the robustness of our fit to the number of red- 
shift slices analysed. This test was performed on the conservative 
data sample, excluding the early SDSS data and the 2dFGRS. In 
our default analysis we fit power spectra calculated for six redshift 
slices, and the resulting likelihood surface for the late SDSS sam- 
ple is shown in panel (b) of Fig. [7] For comparison, panel (e) of 
Fig-H]shows the likelihood surface calculated using power spectra 
from only three redshift slices (details of the slices chosen are pre- 
sented in Section[3](. Because we are only fitting two Dv (z) nodes, 
these should be constrained by our reduced fit using three redshift 
slices. Panel (e) of Fig. [8] shows that this is true, but comparison 
with panel (b) of Fig.|7]shows that the constraints are tighter if we 
model power spectra from six redshift slices. Clearly, extra infor- 
mation is available from the extra redshift slices, and we therefore 
fit to six redshift slices for our default analysis. 

It is interesting to test if there is sufficient information to con- 
strain the shape of Dv(z) beyond our simple spline model with 
2 nodes. Results from fits allowing four Dy(z) nodes are shown 
in Fig. [9] There is a clear maximum in the slices through the like- 
lihood surface close to the ACDM model, but the surface is noisy, 
and there are secondary maxima present. There is a strong degener- 
acy between Dv"(0.3) and _CV(0.45), and between Dv'(O.l) and 
IV (0.2): the data contain limited information to distinguish the 
shape of the distance-redshift relation between these redshifts. Con- 
sequently, we do not try to extract this information, instead concen- 
trating on fits where there are only two nodes in Dy (z). 



8.4 The covariance matrix 

Because we are analysing overlapping shells in redshift, the power 
spectra will be strongly correlated and the estimation of the co- 
variance matrix will be in error if we do not have sufficient mock 
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Figure 7. As Fig. [4] but now considering results from 4 choices of catalogue: (a) all SDSS and non-overlapping 2dFGRS data, (b) excluding both the early 
SDSS data and 2dFGRS, (c) excluding the 2dFGRS, (d) excluding the early SDSS data. 



catalogues. In order to test this, we have recalculated our covari- 
ance matrix using 1/3 as many LN catalogues, and have used this 
matrix to recalculate the required corrections to the confidence in- 
tervals using independent sets of LN catalogues. We find consistent 
results in the factors required to match the confidence intervals to 
those expected for a multi-variate Gaussian distribution. We have 
also performed a full analysis using this reduced covariance matrix, 
and find results consistent with using our default covariance matrix. 



8.5 Calibration 

The likelihood surface shown in panel (f) of Fig. [8] was calculated 
using a SDSS galaxy samp le with luminositi es calibrated using 
the photometric calibration (Tucker et al. 2006), prior to the uber- 
calibration analysis jPadmanabhan et al2008l) . This affects the cal- 
culation of the redshift completeness for any region observed, and 
also the luminosity-dependent weights applied to the SDSS galax- 



ies. The effect of this calibration change on our results is small, 
and there is no significant change between the likelihood surface in 
panel (f) of Fig. [8] and that in panel (b) of Fig. [7] where the uber- 
calibration data set was used. 



8.6 BAO damping scale 

Panel (g) of Fig. [8] shows the likelihood surface if we allow the 
BAO damping scale to be a free parameter in the fit, placing a sim- 
ple Gaussian prior on its value Ddamp = 10 ± 5 h^ 1 Mpc. This 
prior on th e BAO damping scale is conservative. From simulations, 
iReid etal.l f2008) found D damp = 9.2 ± 1 h' 1 Mpc, with no vari- 
ation with redshift for < z < 0.5 for halo density fields, and 
Dda,mp = 9.7±1 h~ l Mpc for density fields matc hed to the LRGs. 
The m ild cosmological dependence suggested by lEisenstein et all 
d2007l) shows that the main cosmological dependence is through 
the linear growth rate; current constraints on as are much better 
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Figure 8. Likelihood contour plots as in Fig. [4] for the SDSS data, excluding the early data and the 2dFGRS, but now calculated for (e) fit to three power 
spectra, (f) old rather than uber-calibration. (g) allowing the BAO damping scale to vary with a simple Gaussian prior Ddamp = 10 ± 5 h Mpc, (h) 
randomized galaxy redshifts used to give the expected radial galaxy distribution. 



than that required to significantly change Ddamp, and we consider 
±5 ft - 1 Mpc to be a conservative prior. Allowing the damping 
scale to vary degrades the constraint, increasing the size of the pa- 
rameter confidence regions. The best-fit solution does not move sig- 
nificantly, suggesting that our default assumption of a fixed damp- 
ing scale is sufficiently accurate to current data precision. 



8.7 Radial galaxy distribution model 

Finally, analysis run (h) shows the constraints if we use a random 
catalogue where we randomly choose a galaxy redshift for each an- 
gular position chosen, i.e. to model the expected redshift distribu- 
tion (n(z)), we sample from the galaxy redshift distribution. This 
test was designed to investigate the dependence of the analysis on 
how well we model the radial galaxy distribution. Randomly sam- 
pling galaxies to obtain this distribution, perfectly matches the red- 



shift distribution of the galaxies and that of the random catalogue 
used to define the survey region. In fact, we see no change in our 
results if we do this rather than using a smooth fit to the redshift dis- 
tribution. This gives us confidence that our results are not sensitive 
to this modelling. 



8.8 Secondary likelihood maxima 

In the likelihood surfaces in Figs. [7] & [8] we see secondary likeli- 
hood maxima, which appear to lie on a degeneracy stretching from 
£>y(0.2) = 700 Mpc, D v (0M) = 1500 Mpc to D v (0.2) = 
600 Mpc, IV (0.35) = 1000 Mpc. These minor peaks in the likeli- 
hood, which appear as isolated islands in the likelihood surface are 
of lower significance than the strong peak close to the parameters of 
a concordance ACDM model. Tests have shown that the secondary 
peaks result from the interplay of two competing effects, which are 
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Figure 9. Contour plots showing slices through the likelihood for four Dv(z) cubic spline nodes at z = 0.1, z = 0.2 z = 0.3 and z = 0.45, calculated 
for our default analysis using six power spectra, uber-calibration, and a fixed BAO damping scale of -Ddamp = 10 h~ x Mpc. Shaded regions are plotted for 
— 21n£/£truc < 2.3, 6.0, 9.3, which for a multi-variate Gaussian distribution with two degrees of freedom correspond to 68%, 95% and 99% confidence 
intervals. Likelihoods were adjusted to match these Gaussian confidence intervals as described in Section|5] In each panel, the nodes that are not shown were 
fixed at the default ACDM (f! m = 0.25, Qa = 0.75) values. We use shaded regions in this plot to show the likelihood surface, compared with the contours 
in Figs.[7]&[H]because the likelihood surface is more complicated with four nodes, and the shading helps to distinguish peaks from troughs. 



themselves a result of using the wrong cosmology to analyse the 
BAO. These are: 



(i) A shift in the BAO position, 

(ii) An increase in the width of the window associated with each 
band-power, caused by BAO in different redshift shells being out- 
of-phase. This can smooth out the BAO signal. 



Secondary maxima are produced where the BAO shift and the 
smoothing "balance". If we redo the analysis ignoring the second 
effect by assuming that the window function is a <5-function centred 
on the peak, these secondary maxima are removed. 
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Figure 10. The expected shift recovered from an analysis of the BAO po- 
sition in a spherically averaged galaxy power spectrum, if there are radial 
and angular distortions induced by assuming an incorrect cosmology when 
analysing the data. The thick solid contour shows no residual shift, while the 
dotted contours show a 1 % shift. For comparison we plot the expected be- 
haviour for an isotropic power spectrum H(z) oc D 2 A , and for an increased 
importance of the radial distortion H (z) oc Da (dashed lines). The top 
panel approximates redshift-space, by weighting the power in the spherical 
average by (1 + fifi 2 ) 2 , with (3 = 0.25, matching that expected for the 
SDSS LRGs, while the bottom panel does not include this weighting. For 
comparison, the thin solid contour in the top panel marks no residual shift 
for data with j3 = 1, showing that we should expect the radial signal to 
increase in importance for such a sample. 



8.9 Dependency on Dy 

A possible concern about our method of analysis is that we as- 
sume a fiducial ACDM model to convert redshifts to comov- 
ing coordinate distances and measure the position of the BAO in 
the spherically-averaged power spectrum. If the true cosmolog- 
ical model has different angular diameter distance-redshift rela- 
tion Da(z) and Hubble parameter H(z) than this fiducial model, 
this would cause angular and radial distortions in the density field 
from which we estimate the power spectrum. By presenting results 
in terms of Dy we remove the anisotropic information, and as- 



sume that the expected BAO position for all cosmological mod- 
els is solely dependent on their predicted value of Dy ■ This must 
break down for models that behave very differently from our fidu- 
cial ACDM model. 

We now test the sensitivity of the assumption that the BAO 
position in the spherically averaged power spectrum only depends 
on Dy for cosmological models that predict significant anisotropic 
distortions in the density field away from our fiducial model. To do 
this, we compute the shifts of the BAO position expected when one 
measures the spherically-averaged power in either real or redshift 
space for such models. To simplify the analysis, we assume that 
the BAO in the spherically averaged P(k) will be shifted by the 
average of the shifts in k predicted over all angles: ie. our BAO fit 
recovers the weighted mean shift in the 3D power. In redshift-space 
we also follow the distant observer approximation, and assume that 
the angular dependence of the true 3D power spectrum is given by 
(1 + /3fi 2 ) 2 , where fi is the cosine of the angle to the line-of-sight 
and P = fi^; 55 jb. The anisotropy in the observed power spectrum 
caused by redshift-space distortions will act as a weight when we 
spherically average. 

For the SDSS LRGs, which provide most of our cosmolog- 
ical signal, we take an effective redshift of z = 0.35, and as- 
sume a ACDM model with Q m (z = 0) = 0.25, giving Q m (z = 
0.35) = 0.45 . The LRGs are strongly biased and the model of 
iTegmark etai] ( 12004) gives an effective relative bias for our sam- 
ple, which we correct for in the power spectrum calculation, of 
(b/b t ) = 1.9. Matching the norm alisation of the measured LRG 
power spectrum (Reidetal. 2009) gives that fe* = 1.34 assum- 
ing the LRG clusteri ng is constant in comoving coordinates (e.g. 
IPercival et ail l2007bh . and that ag (matter, z = 0) = 0.8, so 
cr 8 (matter, z = 0.35) = 0.68 dKomatsu et al]|2009l) . This sug- 
gests that we should expect ft ~ 0.25 for the LRG power spec- 
trum, and we show contours calculated assuming j3 = 0.25 in 
Fig. [10] which we compare with the prediction for f3 — 1. Note that 
our luminosity-dependent weighting means that we are upweight- 
ing highly biased galaxies, and that our analysis will therefore have 
a smaller effective (3 than analy ses without such weighting , such as 
the measurements presented bv lCabre & Gaztana ga (2009). 

Fig. [JJj] shows the relation between radial and angular distor- 
tions, H I Had and Dyi/Dyi.fid, which give rise to zero and ±1% 
shift in the spherical averaged power spectrum. Here, Had is the 
fiducial value of H, and similarly for Da- For general cosmolog- 
ical models, H / Had and Da/ Da,sa will depend on redshift, so 
that the final effective shift will be an average over a trajectory in 
this diagram which is determined by the model to be tested. Fig. 1101 
also shows the expected line of zero average shift we would expect 
if the BAO position only depends on Dy(z), which would lead to 
behaviour such that H(z) oc D\. For comparison, we show the 
prediction for a model with increased importance of the radial dis- 
tortions, with H(z) oc Da- This is included because we would 
expect that the redshift-space distortions will increase the impor- 
tance of the radial information. However, the H(z) oc D\ line is 
a significantly better fit, even in redshift-space. The H(z) oc D 2 A 
line does not cross the contours marking a 1% average shift for our 
redshift-space power spectrum, showing that the assumption that 
the recovered BAO position only depends on Dy at most produces 
a 1% systematic in the best-fit for models with an anisotropy dis- 
tortion away from our fiducial model of up to 20% in the radial 
direction. Such a 1% systematic shift, which requires a model that 
is extremely discrepant from ACDM, is significantly below the sta- 
tistical precision of our 2.7% accuracy distance measurement. It is 
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therefore a reasonable approximation to use our measurements of 
Dy to constrain a wide variety of cosmological models. 



9 COSMOLOGICAL PARAMETER CONSTRAINTS 

We now apply our full constraints to a cosmological parameter 
analysis. We assume that the likelihood of a model is given by 
a multi-variate Gaussian distribution around the Dv (z) measure- 
ments given by Eq. ([3), with the inverse covariance matrix of 
Eq. ((5). Throughout this section we consider four models: a flat 
universe with a cosmological constant (ACDM), a ACDM universe 
with curvature (oACDM), a flat universe with a dark energy compo- 
nent with constant equation of state w (wCDM), and a wCDM uni- 
verse with curvature (o wCDM). This is the same model set consid- 
ered bvlReid et alj d2009h . We use a modified version of COSMOMC 
dLewis & Bridlell2002l) to perform the likelihood calculations. 



9.1 SN + BAO + CMB prior likelihood fits 

We first consider the constraints excluding the angular acoustic 
scale in the CMB, in order to consider data that are independent 
of the dark energy behaviour at the redshifts beyond our sample. 
This is important because it ensures that our results only depend 

on the acceleration of the Universe at late times and so do not 

j ■ 1 1 —\ 

depen d on so-called early dark energy models (Wetterich 1988; 
iRatra & Peebles! 1 19881 ; Iziatev et al.|[l999l ; ISteinhardt et alj|l999h . 
which have non-negligible dark energy at early times. We take 
Gaussian priors Q. c h 2 = 0.1099 ± 0.0063 and Q, b h 2 = 0.02273 ± 
0.00061 from the CMB; these constraints from the ratio of peak 
heights in the WMAP5 data alone do not relax when Q,k and w are 
allowed to vary. We also impose weak priors on —0.3 < fife < 0.3 
and —3 < w < 0. The parame ter constraints from the combi- 
nation of Union supernova (SN) ( Kowalski et al. 2008) and BAO 
likelihoods with these priors are presented in Table [4] The best- 
fit value of Q. m ranges from 0.286 to 0.290, with the 68% con- 
fidence interval, ±0.018, while the mean value of Ho varies be- 
tween 67.8 km s" 1 Mpc -1 and 68.6 km s" 1 Mpc" 1 , and the 68% 
confidence interval remains ±2.2 km s" 1 Mpc -1 throughout the 
four models. In Section [7] we derived BAO only constraints of 
±0.018 on tt m and ±2.2 kms" 1 Mpc^ 1 on H , for fixed Q m h 2 . 
If we include the 4.8% error on n m h 2 from the WMAP5 measure- 
ment, then we should expect these errors to increase to ±0.019 
on fi m and ±2.3 km s -1 Mpc -1 on Hq. These agree perfectly 
with the COSMOMC results if we exclude the supernova data, 
so the small difference between the errors in Tableland those ex- 
pected is caused by the supernova data helping to constrain Q. m and 
Hq slightly. Similarly, the best-fit values of these parameters agree 
for COSMOMC results excluding the supernova data. Comparison 
between Tableland Section[7]shows that the inclusion of the super- 
nova data is moving the best-fit slightly: ±0.004 in Q m and —0.5 in 
H for the ACDM model. The COSMOMC analysis therefore val- 
idates the simple derivation presented in Section]?] In the space of 
models considered here, the BAO constraint on Dv (0.275) already 
restricts 7JV(0.35)/.CV(0.2) to a much smaller region than our 
constraint in Eq. Jilt allows. While the combination of these data 
and our priors are unable to constrain fife, w is constrained at the 
±0.11 level. For the owCDM model, the weak prior on fife leads 
to an apparent constraint on w, but these errors depend strongly on 
the prior. 

The data are compared with the best-fit ACDM model in 
Fig. [TT] Three ways of considering the data constraints are shown 
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Figure 11. The BAO constraints (solid circles with lcr errors), compared 
with the best-fit ACDM model. The three panels show different methods of 
using the data to constrain models. 



in different panels. In the bottom panel we plot Dv(z) / Dv (0-2), 
which corresponds to matching the geometry at z = 0.2 and 
z = 0.35 so the BAO match at these redshifts, without includ- 
ing information about the comoving position of the BAO. In the 
middle panel we plot r s (zd) / Dv {z), where we now have to model 
the comoving sound horizon at the drag epoch. In the top panel 
we include a constraint on the sound horizon projected at the last- 
scattering surface as observed in the CMB. Marginalising over the 
set of flat ACDM models constrained only by the WMAP5 data 
gives r 3 (z d )/S k (zd) = 0.010824±0. 000023, where S k (z d ) is the 
proper distance to the ba ryon-drag redshift Zd = 1020.5, as mea- 
sured by WMAP5 team dKomatsu et~ai1l2009h . Ignoring the negli- 
gible error on this quantity, we combine with the BAO results to 
measure Sk{zd)/ 'Dy(z). This effectively removes the dependence 
on the comoving sound horizon at the drag epoch, anchoring the 
BAO measurements at high redshift: here we have done this at the 
baryon-drag epoch so the CMB constraint has matched sound hori- 
zon and projection distance. 

9.2 CMB + BAO likelihood fits 

We now turn to the constraints from our BAO measurement com- 
bined with the full WMAP5 likelihood, including the constraint on 
r s (zd) /Da at the time of decoupling. While this extra constraint 
can break degeneracies between fi m , fife, and w inherent in our 
BAO constraints, the results are now sensitive to our assumption of 
a constant dark energy equation of state w at z > 0.35. Results for 
the four models are presented in Table[5] 
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parameter 


ACDM 


oACDM 


wCDM 


owCDM 


#0 


0.288 ±0.018 
68.1+H 


0.286 ±0.018 
68.6 ± 2.2 


0.290+g-g" 
— u . u i y 

67.8 ± 2.2 


0.286 ±0.018 
68.2 ± 2.2 


i~l k 
w 

n A 


_ 

0.712 ±0.018 


-0.097 ± 0.081 

n SI 1 +0-084 


- 

-0.97 ±0.11 

710+ 019 


-0 1 99+0 080 

U - 1MM _0.089 

n cqc+0.083 
n O13+0.092 


do. 275 


0.1381 ±0.0034 


0.1367 ± 0.0036 


0.1384 ±0.0037 


0.1386 ± 0.0037 


D^(0.275) 


1111 ± 31 


1120 ±33 


1109 ± 32 


1108^33 


/ 


1.662 ±0.004 


1.675 ±0.011 


1.659 ±0.011 


1.665 ±0.011 


Age (Gyr) 


Uralsk 


14.43 ± 0.48 


13.95 ± 0.36 


14.38 ± 0.44 



Table 4. Marginalized one-dimensional constraints (68%) for BAO+SN for flat ACDM, ACDM with curvature (oACDM), flat wCDM (wCDM), and wCDM 
with curvature (owCDM). The non-standard cosmological parameters are do. 275 = T's(zd)/ Dv (0.275) and / = £V(0.35)/ Dy (0.2). We have assumed 
priors of Q, c h 2 = 0. 1099 ± 0.0063 and n b h 2 = 0.02273 ± 0.00061, consistent with WMAP5-only fits to all of the models considered here. We also impose 
weak flat priors of —0.3 < Q k < 0.3 and —3 < w < 0. 



parameter 


ACDM 


oACDM 


wCDM 


owCDM 


owCDM+SN 


owCDM±H 


owCDM+SN±#o 




0.278 ±0.018 
70.1 ± 1.5 


0.283 ± 0.019 
68.3+2 i 


0.283 ± 0.026 
69.3 ± 3.9 


240+ ' 044 
u.ziv_ 043 

75.3 ±7.1 


0.290 ± 0.019 
67.6 ± 2.2 


240+ - 025 

V.Z4U_ Q24 

74.8 ± 3.6 


0.279 ± 0.016 
69.5 ± 2.0 






-0 007+ 006 




-0.013 ± 0.007 


-0.006 ± 0.008 


-0.014 ± 0.007 


-0.003 ± 0.007 


w 




-0.97 ±0.17 


i cq+0.51 


-0.97 ± 0.10 


-1 49+ 32 
-0.31 


-1.00 ±0.10 


n A 


0.722 ±0.018 


0.724 ± 0.019 


0.717 ±0.026 


0.772 ± 0.048 


0.716 ±0.019 


0.773 ± 0.029 


0.724 ± 0.018 


lOOf^/i 2 

T 


2.267 ± 0.058 
0.086 ± 0.016 


2.269 ± 0.060 
0.089 ± 0.017 


2.275 ± 0.061 
0.087 ± 0.017 


9 9C/1+0.062 
z ' ZJ *-0 .061 

0.088 ± 0.017 


2.271 ±0.061 
0.089 ±0.017 


2 254+ - 061 
Z,ZJ -0.062 

0.088 ±0.017 


2.284 ± 0.061 
089+ - 017 


n s 


0.961 ± 0.013 


0.963 ± 0.014 


0.963 ± 0.015 


0.958 ± 0.014 


0.963 ± 0.014 


0.957 ± 0.014 


0.964 ± 0.014 


ln(10 10 A 05 ) 


o 074 +0.040 
o.uit_ Q39 


3.060 ± 0.042 


3.070 ± 0.041 


3 062+ ' 042 
o.uoz_ Q43 


3 062+ 041 
o.uoz_ Q42 


3.062 ± 0.042 


3.072 ± 0.042 


do. 275 


0.1411 ±0.0030 


0.1387 ±0.0036 


i 4 n4+ 0036 

U.14U4_q 0035 


0.1382 ±0.0037 


0.1379 ±0.0036 


n 1 oov+0.0036 
u.iooi_ 0037 


1402+ ' 0033 


(0.275) 


1080 ± 18 


1110+ 32 


1089 ±31 


1111 ± 33 


1115 ± 32 


1107 ±31 


10911 27 


/ 


1.6645 ± 0.0043 


1.6643 ±0.0045 


1.661 ± 0.019 


1.72 ±0.056 


1.660 ±0.011 


1 71O7+0-0337 
' 10 ' -0.0334 


1.6645 ±0.0107 


Age (Gyr) 
tt c h 2 

^tot 


13.73 ±0.12 
0.1139 ±0.0041 


14.08 ± 0.33 

1090+ 0060 
u.iuju_ 0061 

1 007+ 006 

1 - uu ' -0.007 


13.76±°-i| 

n 1199 +0.0068 
u.lizz_ 00g9 


14.49 ± 0.52 

no7+ 0063 
u.nu(_ 00g2 

1.013 ±0.007 


14.04 ± 0.36 

1096+ 0061 
u.iuao_ 0062 

1.006 ±0.008 


14.48 ± 0.48 
0.1108_ 0061 
1.014 ±0.007 


13.86i°; 34 
0.1115 ±0.0061 
1.003 ±0.007 


&8 


0.813 ±0.028 


0.787 ± 0.037 


n 7q2 +0.081 
U. <y^_ .082 


0.907 ±0.117 


U.78U_ 053 


0.904 ± 0.074 


n am +0.053 
0.80l_ 052 



Table 5. Marginalized one-dimensional constraints (68%) for WMAP5+BAO for flat ACDM, ACDM with curvature (oACDM), flat wCDM (wCDM), 
wCDM with curvature (owCDM), and owCDM including constraints from supernovae. The non-standard cosmological parameters constrained by the BAO 
measurements are do.275 = T s (z t i)/ Dy (0.275) and / = Dy(0. 35)/Dy (0.2). 



For the ACDM model, we find Q m = 0.278 ± 0.018 and 
H = 70.lil.5kms" 1 Mpc" 1 , with errors significantly reduced 
compared to the WMAP5 alone analysis (Q, m = 0.258 ± 0.03 and 
Ho = 7 0.51 2 . 7 kms" 1 Mp c -1 ). Similar limits on Q, m were ob- 
tained by Roz o et all (2009) who used the maxBCG cluster abun- 
dance and weak-lensing mass measurements to similarly break the 
tight WMAP5 constraint on Q m h 2 . 

Fig.[l2]shows the impact of relaxing the flat, ACDM assump- 
tion. The WMAP5 results alone tightly constrain Q m h 2 in all of 
these models (dashed lines), but low redshift information is neces- 
sary to constrain Q m and Ho separately. Allowing w 7^ — 1 relaxes 
the constraint on fl m from the BAO measurement, and in addition 
allowing fl^ 7^ relaxes the constraint even further. The impact 
on the constraints on Q m and Ho is shown in the lower right panel. 
All of the contours lie along the banana with Qmh 2 fixed from the 
CMB. 

In the oACDM model, the combination of scales measured 
by the CMB and the BAO tightly constrain the curvature of the 
universe: fit = — 0.0071o'oq 7 . The constraints on f2 ro and Ho in 
this model are well described by Eqns. jl6\ & H7\ , while in the 



wCDM cosmology they degrade because w is not well-constrained 
by the low redshift BAO information alone. 

When the parameter space is opened to both curvature and 
w, the WMAP5 data are not able to eliminate the degeneracy be- 
tween Q m and w in the BAO constraint. The constraints relax 
to n m = 0.240l^ 4 3 and H = 75.3 ± 7.1 kms" 1 Mpc" 1 ; 
Qk = —0.013 ± 0.007 is still well-constrained but w is not (see 
Fig.ll2t. Including the constraints from the Union Supernova Sam- 
ple breaks the remaining degeneracy, and we recover the tight con- 
straints on Q m and Ho given in Eqns. l |16t & (T7J. These con- 
straints, and the relative degeneracies induced and broken by dif- 
ferent data sets, are shown in Fig. [13] For each of the four models 
considered, the central values for Q. m and Ho change only slightly 
when the full WMAP5 likelihoods are used (Table|5j instead of pri- 
ors on Qth 2 and Q c h 2 in combination with the Union SN sample 
(TableH. 

Table|5]also lists the best fit cosmological age (i.e. time since 
the Big Bang) for different cosmologies and data sets. While the 
age is very well determined for ACDM and wCDM, there is a de- 
generacy between age and curvature that increases the uncertainties 
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Figure 12. WMAP5+BAO constraints on Cl m h 2 , f2 m , and Ho for ACDM (solid black contours), oACDM (shaded green contours), wCDM (shaded red 
contours), and owCDM (shaded blue contours) models. Throughout, the solid contours show WMAP5+LRG ACDM constraints. The first three panels show 
WMAP5 only constraints (dashed contours) and WMAP5+BAO constraints (colored contours) in the Q m h 2 - Q m plane as the model is varied. In the lower 
right we show all constraints from WMAP5+BAO for all four models in the Q m — H plane, which lie within the tight Q m h 2 Ri 0.133±0.006 WMAP5-only 
constraints. 



and allows for an older age in oACDM and owCDM information. 
Adding SN and Ho measurements reduces these uncertainties and 
implies a best fit age of 13.1 



StSJS Gyr. 



9.3 Comparison with lRiess etafl fe009l) H 

iRiess et alj j2009h recently released a new determination of the 
Hubble constant using a differential distance ladder: Ho — 74.2 ± 
3.6 kms" 1 Mpc -1 . This value, as well as the values Ho ~ 
68kms _1 Mpc -1 determined in Table 0] using BAO, SN, and a 
WMAP5 prior on Q, c h 2 and Qbh 2 , are within ~ la of the mean 
value determined from WMAP5+BAO in a ACDM model, 70.1. 
In the wCDM model, combining this new Ho with the WMAP5 
likelihood constrains w = —1.12 ± 0.12. In Table [5] we show 
MCMC results for the owCDM model for WMAP5+BAO+#o and 
WMAP5+BAO+i?o+SlsQ. In this model, the supernova data are 
more effective than Ho at breaking the long degeneracy in the 



1 We account for the small cosmology dependence in the H p constraint 
(seen as a slight degeneracy between Ho and w in fig. 14 of [Riess et alj 
120091) by considering it as a constraint on the inverse luminosity distance at 
the effective redshift z = 0.04 (Riess private comm.) 



WMAP5+BAO constraints. Combining WMAP5+BAO+SN+// , 
the mean parameters are quite close to ACDM: Qfc = —0.003 ± 
0.007 and w = -1.00 ± 0.10, and fi m = 0.279 ± 0.016 and 
Ho = 69.5 ± 2.0kms _1 Mpc" 1 are also well-constrained. 



10 COMPARISON WITH DR5 ANALYSES 

In lPercival et al. I d2007ch . we presented BAO measurements calcu- 
lated from fitting power spectra calculated for three samples drawn 
from the combined SDSS+2dFGRS catalogue, using the SDSS 
DR5 data. The full catalogue was split into galaxy populations, 
rather than redshift slices, corresponding to the SDSS LRGs, the 
2dFGRS+SDSS main galaxies, and the combined sample. From 
this, we obtained the distance constraints r s (zd)/Dv(Q-'2) = 
0.1980 ± 0.0058 and r 3 {z d )/D v (0.35) = 0.1094 ± 0.0033 with 
correlation coefficient 0.39, which gives a distance ratio measure- 
ment of ZV(0.35)/Dv(0.2) = 1.812 ± 0.062. The concordance 
ACDM value is IV(0.35)/IV(0.2) = 1.67, measured using the 
SNLS supernova data, which is discrepant with the published DR5 
BAO results at the 2.4a level. The analysis of mock catalogues pre- 
sented in Section [5] showed that the cubic spline x BAO method 
underestimates the true distribution of recovered distances, given 
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Figure 13. For the owCDM model we compare the constraints from 
WMAP5+BAO (blue contours), WMAP5+SN (green contours), and 
WMAP5+BAO+SN (red contours). Dashed and solid contours highlight 
the 68% confidence intervals for the WMAP5+BAO and WMAP5+SN 
models respectively. 



noisy data, which produce a non-Gaussian likelihood surface. We 
should therefore incre ase the errors on the DR5 measurements of 
IPercival et all J2007ch by at least a factor of 1.14, which is the 
correction derived from the fits to three DR7 power spectra. If 
we do this, the revised DR5 constraints are r s (zd)/Dv (0.2) = 
0.1981 ± 0.0071 and r s (z d )/D v (0.35) = 0.1094 ± 0.0040 with 
correlation coefficient 0.38, which gives a distance ratio measure- 
ment of 7JV(0.35)/LV(0.2) = 1.813 ± 0.073. The discrepancy 
between the old DR5 constraints and the SNLS ACDM value is re- 
duced to < 2<r. Because the DR5 data were noisier than the DR7 
data, we should expect the likelihood surface to be less like a Gaus- 
sian prediction, and the correction actually should be slightly larger 
than that for the DR7 data. 

Of all the changes implemented between this DR7 analysis 
and the analysis of the DR5 data, it was the increase in the num- 
ber of random points used to quantify the survey geometry that 



had the most effect when comparing different catalogues. We now 
find consistent results, given in Table [3] for all catalogues and 
analysis variations presented in Section [8] When translated into 
constraints on the distance ratio, for the full catalogue we find 
£V(0.35)/LV(0.2) = 1.736±0.065. Using only 3 redshift slices 
we find Dv(0.35)/Dv(0.2) = 1.765 ± 0.079. If the 0.5cr differ- 
ence is not due to chance, the difference between these measure- 
ments could be caused by residual non-Gaussian scatter in the band 
powers. A scenario in which this is reduced by including fits to 
more redshift bins would then explain the observed trend. Exclud- 
ing the 2dFGRS and early SDSS data, the constraint is reduced to 
7Jv(0.35)/Dv(0.2) = 1.747 ± 0.070, which is consistent with 
the t ighter constraint using all of the data. 

I Sanchez et alj ( 120091) , who analysed the SDSS DR6 sam- 
ple, speculated that the discrepancy could be caused by the 
IPercival et al. (2007c) analysis fixing the BAO damping scale. 
However, in our current analysis, if we allow the BAO 
damping scale Ddamp to vary, the derived constraints on 
7JV(0.35)/Dv(0.2) does not change significantly from that re- 
covered in our default analysis. The mild discrepancy with ACDM 
does not appear to be caused by fixing the damping scale. The 
change from photometric calibration to uber-calibration has a rel- 
atively minor effect on the distance ratio, which increases to 
Dv(0.35)/LV(0.2) = 1.748 ± 0.074. Fig. ©shows that the ef- 
fect on the BAO of redshift-space distortions caused by the thermal 
motion of galaxies in clusters is similarly small. Linear redshift- 
space distortions propagate the apparent position of galaxies along 
their velocity vector in a way that simply makes the field look more 
evolved than it is; they do not alter the positions of the BAO. 

In conclusion, the significance of the discrepancy with flat 
ACDM models is reduced because of 

(i) analysis of the non-Gaussian nature of the likelihood surface, 

(ii) analysis of more redshift slices, 

(iii) more accurate determination of the galaxy redshift distribu- 
tion. 



11 DISCUSSION 

In this paper we have measured and analysed BAO from the SDSS 
DR7 sample, which represents the final data set observed using 
the original SDSS spectroscopic target selection al gorithm. We 
have fu rther developed the analysis method used by IPercival et all 
(2007c) to analyse the DR5 sample, including a faster method for 
the calculation of the window function (see Appendix [A), linking 
the cosmological model to be tested with the power spectrum band- 
powers measured. This has enabled us to analyse power spectra 
calculated for six rather than three redshift slices, which would not 
have been possible using the old method. 

In Section [6] we have shown how the distance-redshift con- 
straints at z = 0.2 and z = 0.35 can be decomposed into a single 
distance constraint at * = 0.275, and a "gradient" around this pivot 
given by Dy (0.35) / Dv (0.2). This allows us to easily test the con- 
sistency of the ACDM model without having to compare with ad- 
ditional data. For the best-fit flat ACDM model that matches our 
constraint do. 275 = 0.1390 ± 0.0037, we find that our distance- 
ratio measurement of D v (0.35) / D v (0.2) = 1.736 ± 0.065 is 
consistent at the l.lcr level. 

Now that the SDSS-II sample is complete, the importance of 
including the 2dFGRS data is reduced, and the inclusion only de- 
creases the low redshift z = 0.2 distance error by 4%. As we 
showed in Section [8T2l the inclusion of the 2dFGRS galaxies does 
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not lead to the discrepancy with the ACDM model: including the 
2dFGRS brings our constraint slightly more into line with the pre- 
dictions of ACDM models. 

Of the cosmological parameter constraints presented in Ta- 
bles [4] & [5] perhaps the most impressive are the constraints on 
Q m and Ho. For ACDM models, fitting to BAO and Supernovae 
with priors on Q, m h 2 and Qth 2 gives Ho to 3.2% and fi m to 6.4%. 
These constraints are robust to the behaviour of the Universe at 
high redshift, as they are based only on the distance-redshift rela- 
tion at redshift z < 0.35: we can allow fit 7^ and w 7^ — 1 with 
minimal effect. This weak dependence on w and Qk was shown in 
Eqns. (16) & (TTJ for the BAO data. 

If we allow for the flatness constraint to be relaxed, then 
we obtain Qk ~ —0.007 ± 0.007 from the combination of 
BAO+WMAP5 data. A tight constraint was similarly obtained on 
w = —0.97 ± 0.17 if we relax the A constraint. If we allow both 
the curvature and the dark energy equation of state to vary, we must 
include more data to continue to break the degeneracy between the 
two parameters. We do so by including results from the Union SN 
dataset, giving us D. k = -0.006 ± 0.008 and w = -0.97 ± 0.10, 
consistent with a flat ACDM model. If one allows only w 7^ — 1 OR 
Qk 7^ 0, then the combination of CMB, supernova and BAO data 
has an internal cross-check: opening two degrees of freedom from 
flat ACDM yields results that are consistent with flat ACDM. We 
have also shown that our c onstraints are con sistent with the recent 
re-determination of Ho bv lRiessetai] f2009), and that combining 
this constraint with WMAP5, BAO, and SN in a model where both 
curvature and w vary yields mean parameter values very close to 
ACDM. 

In a companion paper dReid et alj2 009), we consider the LRG 
sample in more detail. The LRGs are distributed in haloes in a 
simple way and we are able to extract the halo power spectrum 
from the data. In addition to fitting the BAO in this power spec- 
trum, we are able to extract limited information about the shape 
of the power, which gives complementary constraints. A detailed 
comparison between the results from our fit to the BAO in red- 
shift slices, performed in a cosmology model-independent way and 
including l ow-re dshift galaxies, with the halo power spectrum of 
iReidetal .1 (2009) is presented in that paper, where excellent agree- 
ment is demonstrated. The data sets are correlated so they should 
not be used together to constrain cosmological models. 

Our analysis highlights the importance of BAO as a key 
method for investigating cosmic acceleration, and shows that the 
method can already provide interesting cosmological constraints. 
Ongoing spectroscopic surveys aiming to use BAO to analyse 
dark energ y include the Baryon Oscillation Spectroscopic Sur- 
vay (BOSS: ISchiegel et al.l2009ah. the Hobby-Eberly Dark Energy 
Experiment (HETDEX; Iffill et al.|[2008h and the WiggleZ survey 
dGlazebrook et alj|2007t) . There are also plans for future surveys 
covering significantly larger volumes of the Universe, and therefore 
observing the BAO signal with higher precision such as the Square 
Kilometer Array (SKA: www.skatelescope.org), and the 
Joint Dark Energy Mission (JDEM: jdem.gsfc.nasa.gov) 
and European Space Agency Euclid satellite mission concepts, 
or the Big Baryon O scillation Spectroscopic Survay (BigBOSS; 
ISchlegel et al] |2009b). Photometric surveys such as the Dark 
Energy Survey (DES: www.darkenergysurvey.org), the 
Panoramic Survey Telescope & Rapid Response System (Pan- 
Starrs: pan-starrs.ifa.hawaii.edu) and the Large Syn- 
optic Survey Telescope (LSST: www . lsst . org) will find BAO 
using photometric redshifts. All of these surveys will measure BAO 
at higher redshifts than those analysed in our paper using SDSS-II 



data: if dark energy does not have a simple explanation, then com- 
parison between future high redshift results and our current under- 
standing of the low-redshift Universe from SDSS-II will provide an 
interesting test of these models. 
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APPENDIX A: CALCULATION OF THE WINDOW 
FUNCTION 

In this Appendix, we describe the method used to calculate the 
mapping between the power spectra in the "true" cosmology to 
be tested, and the measured, or observed, power spectra where 
a ACDM model was used to convert redshifts to distances. This 
window function includes both the effect of the survey geometry 
and the mapp i ng betw een cosmological models. As described by 
iPercival et al. (2007c), we should expect the observed power spec- 
trum to be a convolution of the true power spectrum with a window 
function. 

P{k) ohs = J dk'W(k,k')P(k') truc . (Al) 

The goal of this section is to introduce a fast method by which 
W(k, k') can be calculated for any model. 
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Figure Al. The redshift dependence of galaxy pair-weights for the SDSS 
DR7 LRG and main galaxy samples, and from the combination of the 
two. These curves were calculated assuming a flat ACDM cosmology with 
tt m = 0.25, h = 0.72, & U b h 2 = 0.0223. 



In lPercival etal I d2007d) . this window function was calculated 
using Monte-Carlo realisations of Gaussian density fields, created 
assuming the cosmological model to be tested. These fields were 
then distorted as if they had been analysed assuming a ACDM 
model, and the power spectrum was calculated and compared with 
that input. Using a large number of simple input power spectra, 
we were able to construct the window function from this compar- 
ison. This procedure required significant computational resources 
as many density fields were needed in order to accurately measure 
the window function, limiting the number of models that could be 
tested. In particular, we were only able to consider cubic spline 
models of Dy{z) with two nodes to three power spectra. With a 
faster window function calculation, we can include more nodes, 
and fit to more power spectra. 

For a survey covering a thin shell, the window function relat- 
ing true and observed power is an offset delta function 



W(k,k') = S D [k/k' - e], 



(A2) 



where e = d p (true)/d p (obs) is the ratio of proper distances in the 
true and observed cosmologies. Here we are simply stretching the 
true survey prior to measuring the power spectrum. 

The obvious extension to surveys over a range of redshifts is 
to split the sample into i redshift shells, and to approximate the 
window function as 



W(k,k') = ^6 D [k/k' - u]wi, 



(A3) 



where Wi is the weighted number of galaxy pairs in redshift shell i. 
Because we are now considering a broad survey, this pair weight 
is a function of pair separation. In this paper, we bin pairs of 
galaxies with comoving separation 90/i~ 1 Mpc < (^acdm < 
130 /i" 1 Mpc, where ^acdm is the comoving distance in the 
ACDM cosmology used to convert galaxy redshifts to distances. 
The bin size was chosen to approximately match the BAO scale. 
For the SDSS LRG, main galaxy and combined samples, the galaxy 
pair-weights are shown in Fig. IA1I We also need to allow for dif- 
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Figure A2. Galaxy pair-weights for the SDSS DR7 LRG and main galaxy 
samples, and from the combination of the two, as a function of comoving 
distance shifts. These were calculated assuming that a flat ACDM cosmol- 
ogy with n m = 0.25, h = 0.72, & fl b h 2 = 0.0223 was used to anal- 
yse the data, while the BAO are present in a true cosmological model with 
distance-redshift relation defined by a cubic spline in Dy(z) with nodes at 
z = 0.2 and z = 0.35, with amplitude as shown in the plot. 



ferences in the orientation of galaxy pairs, as the distribution of 
6i should allow the galaxy pairs to be of all orientations. Includ- 
ing radial separations introduces an asymmetric convolution for e;, 
and we have found that this needs to be included in order to pro- 
vide approximately the correct window function shapes. Note that 
Eq. dA3t is exact when there is a perfect dilation of scale between 
the true and observed cosmologies: such stretching of the windows 
can be perfectly represented by this equation. 

For each "true" cosmology to be tested, we can calculate the 
shift in scale that stretches each pair of galaxies because we do 
not measure BAO using this model. We have to allow for the 
angular shift caused by a change in Da{z) and the radial shift 
caused by the true and observed H(z) being different. An exam- 
ple of the weighted distribution of "shifts" expected for a model 
cosmology defined by a cubic spline in Dy{z) with two nodes 
at z = 0.2 and z — 0.35 is shown in Fig. IA2I Here the true 
cosmology has a distance-redshift relation given by a spline fit to 
D v {z), with nodes D v (z = 0.2) = 550 h' 1 Mpc, and D v {z = 
0.35) = 1080 h' 1 Mpc. The ACDM values are D v {z = 0.2) = 
568 h' 1 Mpc, and D v (z = 0.35) = 949 h' 1 Mpc, so at redshift 
z = 0.2, BAO in the true cosmology are stretched to larger scales 
by the analysis method, while those at redshift z — 0.35 are com- 
pressed to smaller scales. For the SDSS main galaxies, with median 
redshift close to z ~ 0.2, dtrue/rfobs < 1, while for the LRGs, with 
median redshift z ~ 0.35, dtrue/dobs > !■ 

For each "true" cosmological model, the window function re- 
lating the true and observed power spectra was calculated by con- 
volving the standard window function for the ACDM model, by 
the distribution of shifts such as that shown in Fig. IA2I For the 
models shown in Fig. lA2l we have calculated the window function 
using the approximate method outlined in this Appendix , and us - 
ing the Monte-Carlo method described by IPercival et alj j2007ch . 
A comparison of the windows is presented in Fig. IA3I Reasonable 
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Figure A3. Window functions for three values of k, calculated for the SDSS LRG, main galaxy and combined catalogues. Dotted lines represent the windows 
for our fiducial ACDM cosmology. Solid and dashed lines show the window functions, if the true cosmology were different, but the data were analysed 
assuming that the fiducial ACDM cosmology i s correct. The solid l ines were calculated using the procedure outlined in this Appendix. Dashed lines were 
calculated using the Monte-Carlo procedure of lPercival eTaT] fe007ch . 



agreement is found between the different methods: it is clear that 
the approximate method of splitting into shells recovers the main 
features of the window function. The agreement is not perfect, as 
expected given the approximate nature of our calculation. Because 
we analyse the data using a ACDM model, the window will be 
correct for this model, and will only deviate if we consider signifi- 
cantly different distance-redshift relations. 
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